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Abstract 

Extensions of a previously presented Landau-de Gennes type liquid 
crystalline phase transition model for the direct isotropic/smectic-A (lamel- 
lar) a phase transition to the direct isotropic/smectic-C (tilted lamellar) 
transition are studied. Two different proposed extensions to the model 
are studied both in the context of ideal ordering (point volume) and full 
three-dimensional uniaxial scalar/vector decomposition. Recommenda- 
tions based upon the inclusion of essential physics and computational 
feasibility are made, distinguishing each of the two proposed extensions 
based upon these criteria. Additionally, it is found that the approach used 
for the isotropic/smectic-A model to deterministically compute phase dia- 
gram data is not possible for either of the two proposed isotropic/smectic- 
C models. 
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1 Introduction 

Liquid crystalline phase ordering or self-organization is a pervasive phenomena 
observed throughout the physical world. The study of this class of soft matter is 
becoming increasingly more prevalent due to liquid crystal-based technological 
advances in the past few decades. Apart from the development of new types of 
display technology and high-performance materials, much interest is focused on 
mesophases (liquid crystalline phases) in Nature. In order to accomplish this 
task, fundemental understanding of mesophases must be extended beyond the 
simplest types of liquid crystals, nematics, where some degree of orientational 
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order is present. Higher order smectic and columnar mesophases, where addi- 
tionally there is some degree of translational order, have been found in a vast 
number of biological materials and processes. Harnessing the self-organizing 
properties of liquid crystals to solve engineering problems can thus be extended 
from the current state of the art, manipulation of bulk material properties, to 
more complex techniques such as biosensing and biomimetics . 

One of the main challenges inhibiting the study of smectic liquid crystals are 
the time (ns) and length (nm) scales of their molecular organization. As a re- 
sult, experimental approaches are mainly limited to capturing static phenomena. 
However, theoretical approaches have been able to both enhance experimental 
techniques and independently uncover unknown physics. Numerical simulation, 
harnessing recent computational advances, are at the forefront of these theo- 
retical approaches. In this context, simulations based on continuum models 
are the most viable approach to studying multiple scales, both experimentally 
accessible and inaccessible. 

The goal of past and ongoing work studying the most simple of the smec- 
tic phases, the lamellar smectic-A mesophase, has shown great promise [1,2]. 
Applying similar approaches to the study of a more complex smectic, the tilted 
lamellar smcctic-C mesophase, is therefore a clear progression. A first step to 
modeling and simulation of these and other materials is the development of a 
suitable model and the incorporation of existing experimental data into that 
model, when possible. To achieve this goal the problem can be approached 
through the development of a continuum model that has characteristics con- 
ducive to representing the observed physical system. In this work, two different 
extensions of a Landau-de Genncs type model mesophase transition model [3 6] 
will be analyzed from this viewpoint. The objective of this work is to determine 
the suitability of the two models for use in both phase diagram computation and 
full three-dimensional simulation of the isotropic/smectic-C mesophase transi- 
tion. 

The paper is organized in this way: first a brief introduction will address 

relevent liquid crystal mesophases and the theortical approach used to charac- 
terize them. The Landau-de Gennes type phase transition model for the direct 
isotropic/smectic-A model will then be presented and features explained in the 
context of these types of theoretical approaches. The extensions of this model 
to the isotropic/smectic-G transition will then be presented and discussed. Fi- 
nally, conclusions regarding the suitability of each extension will be made based 
upon both the extent to which each model incorporates essential physics and 
feasibility of simulation. 

2 Background 

2.1 Liquid Crystal Mesophases 

Liquid crystals or mesophases are partially ordered materials that flow like liq- 
uids but have some degree of orientational/translation order like crystals. They 
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are composed of anisotropic molecules called mesogens which can be discotic 
(disc- like) or calamatic (rod- like). Thermotropic liquid crystals are typically 
pure-component compounds that exhibit mesophase ordering most greatly in 
response to temperature changes. Lyotropic liquid crystals are mixtures of meso- 
gens, possibly with a solvent, that most greatly exhibit mesophase behavior in 
response to concentration changes. Effects of pressure and flow fields/bulk ef- 
fects also influence mesophase behavior. This work will focus on temperature 
effects on thermotropic calamatic liquid crystals. 

An unordered liquid where, there is neither orientational nor translation or- 
der of the molecules, is referred to as isotropic. Of the many different mesophases 
observed, the three of interest in this work are the nematic, smectic-A, and 
smectic-C. These mesophases have increasing degrees of liquid crystalline or- 
dering. The nematic phase exhibits purely orientational order where molecules 
conform to an average molecular axis, the nematic director. Smectic-A ordering 
has orientational order as with nematics but, in addition, has one-dimensional 
translational order in the direction of the average molecular axis. This can 
be thought of in terms of two-dimension fluid layers stacked upon each other. 
Finally, the smectic-C mesophase is identical to the smectic-A, except the av- 
erage molecular axis forms an angle with the normal to the layers. Schematic 
representations of these four phases are shown in Figure [TJ Theoretical char- 

(a) (b) (c) (d) 

Figure 1: Graphical representation of mesophase ordering: a) isotropic b) ne- 
matic c) smectic-A and d) smcctic-C 

acterization of mesophase order is accomplished using order parameters that 
adequately capture the physics involved. These order parameters typically have 
an amplitude and phase associated with them. In order to characterize the par- 
tial orientational order of the nematic phase, a second order symmetric traceless 
tensor can be used [3]: 

Q = S (^nn - i/j + (mm - II) (1) 

where n/m/1 are the eigenvectors of Q, which characterize the average molecu- 
lar orientational axes, and S/P are scalars which represent the extent to which 
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the molecules conform to the average orientational axes [7-9] . Uniaxial order is 
characterized by S and n, which correspond to the maximum eigenvalue(and its 
corresponding eigenvector) of the Q-tensor, S = ^^n- Biaxial order is character- 
ized by P and m/l, which correspond to the lesser eigenvalues and eigenvectors, 
P = -| - l^i)- 

The smectic-A and smectic-C mesophases have one-dimensional translational 
order in addition to the orientational order found in nematics. Characterizing 
these mesophases can be accomplished through the use of primary (orienta- 
tional) and secondary (translational) order parameters together [10]. This is 
accomplished using the tensor order parameter ([T]) and the complex order pa- 
rameter [3]: 

^ = ^pe'* (2) 

where (j) is the phase, ip is the scalar amplitude of the density modulation. The 
density wave vector, which describes the average orientation of the smectic- 
A/smectic-C density modulation, is defined as a = V(^/|V0|. The smectic scalar 
order parameter characterizes the magnitude of the density modulation, and 
is used in a dimensionlcss form in this work. In the smectic-A mesophase the 
nematic director and the density wave vector prefer parallel orientation. In the 
smectic-C mesophase a tilt angle 9 is observed between the nematic director 
and smectic wave vector. 

The isotropic/smectic mesophase transition is first-order, where the materi- 
als exhibit discontinuities in their physical properties across the phase transi- 
tion. These transitions have a latent heat associated with them and a regime 
in which both the disordered and ordered phases coexist in a metastable state. 
This phase transition mechanism involves nucleation and growth. Second order 
transitions are have a continuous evolution of the material properties and the 
growth mechanism is through spinodal decomposition. 

2.2 Landau-de Gennes type model for the direct isotropic/ 
smectic-A transition 

A Landau-type free energy expansion is a Taylor series expansion of the free 
energy of the system. The terms included in this expansion are have a great 
effect on both the phenomena that can be captured and the computational 
difficulty in applying the model. These two competing goals along with the 
use of multiple non-scalar order parameters complicates the matter. In this 
particular model three types of terms can be distinguished from each other: 
homogeneous, gradient, and coupling terms. Only terms up to the fourth order 
are included [11] and in addition to this it is assumed that in the vicinity of the 
transition only the coefficients of the second order terms vary significantly with 
temperature [12]. 

A Landau-de Gennes type model for the first order isotropic/smectic- A phase 
transition is used based upon a model initially presented by Mukherjee, Pleiner, 
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and Brand [3, 11]: 

/5-/0 = la{Q:Q)-^b{Q-Q):Q + ^c{Q:Qf + ^a\^\^ + ^(3m* 
-ijV' (Q : Q) - : (Vvl/) (V**) 

+ ^'1 (Vg)' + \l2 (V • Qf + ^hQ : (VQ : VQ) 
+ ^&i|V*|' + ife2|v2^'|' (3) 



A = aQ{T -Tni) 
a = UQiT-TAi) 

where / is the free energy density, /o is the free energy density of the isotropic 
phase, terms 1-5 are the bulk contributions to the free energy, terms 6-7 are 
coupUngs of nematic and smectic order; both the bulk order and coupling of 
the nematic director and smectic density-wave vector, respectively. Terms 8- 
10 are the nematic and terms 11-12 smectic elastic contributions to the free 
energy, respectively. The order parameters are defined in (IT][2|), T is temper- 
ature, Tni/Tai are the hypothetical second order transition temperatures for 
isotropic/nematic and isotropic/smectic-A mesophase transitions, and the re- 
maining constants are phenomenological parameters. Further explanation and 
justification for the use of this high-order model can be found in [1]. 

The first extension to the isotropic/smectic-A model eqn. |3]to the isotropic/smectic- 
C model was proposed by Mukherjee, Pleiner, and Brand [4] and then later 
modified [5]: 

fc-fo = fs + If {Q Q) - VvfV** + h{Q: VVvf) (Q : VVvf*) (4) 

Later, Biscari,Calderer, and Terentjev [6] proposed an alternate extension 
to eqn. O 

fc'fo^fs + If {Q ■■ VvfVvE-*)' (5) 

3 Results and Discussion 

3.1 Homogeneous volume comparison 

Eqns. [4] and [5] can be simplified for use in phase diagram computation by the 
assumption of a homogeneous volume. This homogeneous free energy equation 
is a function of four scalar components of the two original order parameters 
related to the degree of orientational order, positional order, the magnitude of 
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the density wave, and the tilt angle formed between the density wave and the 
nematic director. Thus, following ref. [4], the smectic density wave is fixed along 
the z-axis and the nematic director n takes the form: 

n = cos{6)z + sin{9)y (6) 

where the unit vector z is chosen parallel to the wave vector and y is an arbritary 
orthogonal vector to z. Note that gradients of 5* does not vanish due to the 
inherent gradient of the density wave. The simplified free energy equation for 
the extension proposed by Mukherjee et al [4,5] is: 

/-/o = \aS^ - ^bS^ + \cS' + \a^^ + \f3^' 

1 ■ i2 ..2 4 / „„„2 " 1 



+ \hi^\^ + \b2^\^ (7) 

where g = 27r/(i is the magnitude of the wave vector and d is the smectic layer 
spacing. The minima criteria for the tilt angle 9 involves the partial derivative 
of the homogeneous free energy with respect to 9: 

^ = = eSi^'^q^ cos 9 sin 9-fS'^i;^q^ cos 9 sin 9-hS^i;^q^ ( cos^ 9-- 
oO \ 3 

(8) 

Assuming the smectic-C phase is present S,ip,q,9 > and eqn. [8] simplifies to: 

= e - /5 - hSq"^ (^cos^ ^ ~ 

Clearly the competition of the e//-terms with the /i-term can result in a finite 
value of the tilt angle 9, but the addition of the / term is not higher order in 9 
(compared to the e-term) , thus can be neglected if it does incorporate necessary 
physics (determined in the next section). The simplified free energy equation 
for the extension proposed by Biscari et al [6] is: 

-^S^i;^ - ^eS^\^ ( cos2 9-];]+ IfS^^/q^ ( cos^ - ' 

+ h,ij\' + i&sV'g' (10) 
which results in a minima criteria for the tilt angle 9: 



= eS'^^9^cos6lsin6'-/S'^V' 9 (cos^e*-- )cos6isin6l (11) 
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Assuming the smectic-C phase is present, eqn. [TT]can be simpHfied to: 







fSip'-'q'' ( cos^ e-l 



(12) 



The competition of the e- and /-terms in this model also results in a finite value 
of the tilt angle 9. 

3.2 Scalar/vector decomposition 

Both models can now be decomposed into scalar/ vector form, using the assump- 
tion of uniaxial orientation such that the order parameters are S the uniaxial 
nematic order parameter, n the nematic director, tp the smectic scalar order pa- 
rameter, and a = V(f> the smectic wave vector. Focusing on the nematic/smectic 
coupling terms e, /, and h, the model of Mukherjee et al [4,5] results in: 



-fs^ 

6^ 



nn : (V^Z-VV- + aa) - ^ ((Vi/-)^ + ^^{a^)) 
nn : 

o 



nnnn :: F + -SS T — - {nnS + 5nn) :: F 
9 3 



r = WipVVip + ■i{aV^p){ay^P) + (VV^VV^) 



(13) 
(14) 



clearly both the e- and /-terms are first-order projections of the nematic di- 
rector onto the smectic wave vector. The inclusion of the e term includes the 
essential physics of the first-order projection, and thus the / term can be ne- 
glected. The higher-order h term is necessary in order to have a finite tilt angle, 
but the physical significance of the higher-order director/ wave- vector projection 
is not clearly essential to the model. This term includes coupling of the nematic 
director to the layer curvature, in addition to accounting for both average and 
Gaussian curvature of the smectic layers. Average curvature of the layers is 
accounted for in the existing 62 term (eqn. [3]) and the inclusion of Gaussian 
curvature would have a negligible on model predictions. Additionally, the phe- 
nomenological coefficient itself would be infeasible to measure experimentally 
or to determine through molecular simulation. The model of Biscari et al [6] 
results in a simpler set of nematic/smectic coupling terms: 
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s5 



nn : 

O 



nn : (VV-VtA + aa) + - ((VVj)2 + ^^(a^)) 
o 



(15) 



where now the /-term is simply second-order in the e-term. This simple ap- 
proach both allows for a finite tilt angle and ease of numerical simulation. 
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4 Conclusion 



A theoretical study of two previously proposed models for the direct isotropic/smectic- 
C phase transition [4 6] were studied in the context of ideal ordering (point vol- 
ume) and full tlirec-diniensional uniaxial scalar/vector decomposition. Based 
upon the inclusion of essential physics and computational feasibility, the model 
of Biscari et al [6] was found to be the most computationally feasible of the 
two models. Nonetheless, past approaches to deterministically compute phase 
diagrams for these types of models [2, 13] would not be possible due to the 
additional couplings introduced through the tilt angle. 
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